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ABSTRACT 


An existing computer code, IDARC, is enhanced to permit the modeling of steel moment- 
resisting frames (SMRFs) with the potential for weld failures at beam-to-column connections. 
The steel member model is derived from flexibility formulations in order to allow complex 
degrading hysteresis behavior to be incorporated. A panel zone element is developed to account 
for inelastic shear deformations in the beam-to-column connection region. Finally, a new 
conceptual hysteresis model is developed to represent the force-deformation characteristics at a 


welded connection, before and after weld failure. 


The new models were validated using experimental data from available component tests and an 
existing computer program, DRAIN-2DX. The results of the study indicate that the enhanced 
program, referred to as IDASS, is capable of adequately reproducing observed behavior of 
SMRFs and can be used as an effective tool to investigate the effects of weld failure in steel 


structures under earthquake loading. 


Keywords: computer program; earthquakes flexibility formulation; hysteresis, modeling; steel 


frames, weld fracture 
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1 INTRODUCTION 


The ability of welded connections in steel moment-resisting frames (SMRFs) to dissipate 
energy in the post-yield range has come into question following the failure of thousands of these 
connections in the recent Northridge earthquake. Additionally, the fact that the beams or panel 
zones in the region of the weld failures experienced no inelastic behavior is indicative of a serious 


problem that merits further investigation. 


The existence of innumerable SMRFs with identical connections in earthquake-prone 
regions highlights the urgent need for a systematic evaluation of typical steel frame buildings in 
which the welded connection is modeled as accurately as possible. A reliable estimate of the 
margin of safety against failure of the structure is impossible if the connection behavior is not 


modeled adequately. 


Existing computer programs cannot readily model the behavior of welded connections, 
particularly the effect of sudden weld failures at beam-column connections on overall frame 
response. The purpose of this research effort is to develop suitable element and material models 
that permit the analysis of SMRFs, before and after the failure of welded connections, under 
earthquake loading. Member models will be developed from concepts of distributed flexibility 
rather than existing procedures using concentrated plasticity. The modeling will also include 
consideration of the joint panel region that may experience inelastic behavior and contribute 


significantly to the overall interstory deformation of the building. 
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2 MODELING ENHANCEMENTS 


The IDARC (Kunnath et al., 1992) computational platform was used to carry out the 
following modeling tasks to enable detailed inelastic analysis of SMRFs with or without welded 


connections. 


i Develop a new member model for steel elements using flexibility formulations with ability 
to include finite plastic hinge length at yielding sections. 

a Develop a panel model to account for elastic and potential inelastic deformations of the 
beam-column joint region. 

i Develop a new hysteretic model for steel sections in which the potential effects of weld 


failure is incorporated. 


Details of the aforementioned developments are detailed in the the following sections. 


| 2.1 =Flexibility-Based Member Model for Steel Sections 


Existing formulations for nonlinear analysis of steel frame structures is based on the two- 
component model introduced by Clough et al. (1965) in which the member is subdivided into two 
fictitious parallel elements, one elastic-perfectly-plastic and the other elastic. The first accounts 
for yielding while the second introduces strain hardening. The member stiffness matrix in this 
case is simply the sum of the stiffness of the two components. The computer program DRAIN- 
2D uses this approach, and has been selected as the companion program in this study for 


comparing the results of the analysis using the flexibility-based member model. 


The idea of a discrete hinge length to model the spread of plasticity was first proposed by 
Soleimani et al. (1979). Here, a finite length representing the plastified zone is allowed to spread 
into the member from the ends of the element, while the interior segment remains elastic. A 
similar model was later used by Meyer et al. (1983). More recently, a girder "superelement” 


consisting of Soliemani's spread plastic element in which the length of the plastic zone is varied as 


a function of load history, a joint subelement to account for fixed-end rotations at the beam- 
column interface, and an elastic sub-element to characterize the element behavior prior to 


yielding, was developed by Filippou and Issa (1988). 


Figure 2-1 shows a typical member model that forms the basis of developing the element 
flexibility matrix. The segment AB refers to the clear span of the member so that the lengths A’A 
and B'B are rigid end zones corresponding to the beam-to-column joint. The joint rotations are 
distinguished from the local member deformations at the potential plastic hinge zone. While the 
solution of the structure stiffness matrix is established at the degrees-of-freedom corresponding to 
the center of the joint, the inelastic behavior of each member is a function of the local rotation 
within the clear spans. In Figure 2-1, a member in double curvature is shown. It should also be 
pointed out that the element formulation is not influenced by the incorporation of inelastic joint 
panel deformations since the member model is derived for the clear span length. The introduction 
of joint shear deformations results in an additional degree-of-freedom to distinguish beam and 


column rotations, as described later. 


The incremental moment-rotation relationship is established from the integration of the 


(M/EI) diagram. The flexibility matrix is expressed in the following incremental form: 


AO , iT fir Sa AM, 
lass) “tL fad tow = 


where AO, and A®, are the incremental rotations corresponding to the moment increments 
AM , and aw, , f 4 are the flexibility coefficients, and L, is the clear length of the member. 


Note that the shear deformations can also be directly incorporated into the above formulation. 


Two variations of flexibility are used in the present model, as shown in Figure 2-2. The 
first model considers a linear variation of flexibility. The flexibility coefficients for this case are as 


follows: 
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Figure 2-1. Member Model with Rigid End Zones 


Model 1 


Figure 2-2. Distributed Flexibility Model with and without Spread of Plasticity 


fy = 1/4(EDa + 1/12(EDb (2-2a) 
fo =fy = —-1/12EDa - 1/12(EDb (2-2b) 
fy = 1/12(EDa + 1/4(EDb (2-2c) 


The second model includes the spread of plasticity across a finite hinge length. The 
flexibility distribution in this case assumes the form shown as Model 2 in Figure 2-2. The 


following flexibility coefficients are obtained: 


f,, = B+3a-307 +0° -B-B* -B*> +208 -a’B +087) /12(ED, 


(2-3a) 
+(1-30 +307 -0° +8 +B? +B* —208 +a°B -aB*)/12(El), 
fy =f, =A+at+a’ -—a’? -B-B? +B’ +a°B -oB’)/12(ED), oa 
+(1-a—-a? +0°+BR+B? —B* +aB* —a’B)/12(ED), 
fo =A+a+07 +a° -38 +387 —-B* —208 -a’B +08*)/12(E/), 
(2-3c) 


+(3-a-a7 —a° +38 -3B? +B* +208 +a7B —oB’) /12(EN), 


In the present formulation, the finite hinge length can be specified in two ways: either as a 
fixed quantity expressed as a percentage of the member length or allowed to vary as a function of 
the end moments. In the latter case, the hinge lengths are set to zero during the initial elastic 
phase. Yielding at either end of the member results in a new stiffness matrix, which is then 
constantly updated as the hinge length increases. The hinge length is a function of the previous 
maximum moment and does not change until the "plastic zone" is exceeded by additional inelastic 
excursions. The flexibility matrix for a member needs to be updated for one or both of the 
following reasons: (a) a transition in stiffness as prescribed by the hysteretic force-deformation 


model; and (b) a change in the plastic hinge length. 


22 Joint Panel Model 


The intersection of large cross-sections in the joint region of SMRFs results in a sizable 
panel zone that can deform in shear and contribute significantly to the overall joint rotations. 
Moreover, the shear stresses in the joint may exceed its yield limit leading to hysteretic energy 
being dissipated by the panel. The effect of such inelastic action in the panel zone may alter the 


dynamic response of the overall structural system. 


Existing procedures in frame analysis assumes either (a) the panel is rigid, in which case 
the angle between adjacent members (beams and columns) remains fixed even after the panel zone 
has undergone severe shear deformation, or (b) a linear, elastic relationship exists between the 
shearing forces and panel-zone distortions. In the former case, a single moment and associated 
joint rotation is used at the center of the panel. The latter approach recognizes the significance of 
joint deformations but is incapable of accounting for large inelastic rotations that may occur if the 


yield shear stress of the joint is exceeded. 


A macroelement model to account for additional shear deformations in the joint region is 
developed. The joint region is idealized as a panel zone characterized by pure shear deformations. 
The formulation assumes that the columns above and below the panel, and the beams to the left 
and right of the panel, are connected by rigid links which are capable of independent rotations. 
The resulting formulation adds an extra degree-of-freedom at each node. The inelastic shear- 
deformation characteristics of the joint panel are prescribed by a bilinear hysteretic model based 


on experimentally observed behavior. 


Figure 2-3 shows a typical beam-column joint region with the panel zone. With reference 


to the figure, Mp, 8p, Mc, 9¢ are the moments and rotations of beams and columns, respectively. 


The shear distortion of the panel, y , is the relative change in the rotations of the beam and 


column element, as follows: 


1p =@, -9-) (2-4) 
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(b) Panel Forces and Stresses 


Figure 2-3. Joint Panel Model 


A relationship of the following form can be derived: 


M. =-VG I -I1| (9, 75 
Mile) fe i 0, Se 


where V, = volume of the panel given by (/,b,t, ) where h, is the panel depth, b, is the panel 
width, t, is the panel thickness, and G = shear modulus of the material. The shear vs. shear 
strain behavior is specified by means of a bilinear nondegrading envelope, shown in Figure 2-4. 
The motivation for the loop behavior is derived from observed experimental response of panel 


zone deformation. 


2.3 Hysteretic Model for Steel Sections 


Two separate models were developed for steel sections. The first is a simple bilinear 
model which is commonly used in nonlinear analysis of steel structures. The second model 
attempts to simulate the failure of a welded connection. A brief description of the two models is 


_ provided below. 


2.3.1 Nondegrading Bilinear Model 


The load-deformation path in this model is prescribed by a primary stiffness: k, for the 
loading and unloading segments, and : kz (+ and -) for the strain-hardening or post-yield 
stiffness. The expected load-deformation behavior is shown in Figure 2-4. The yield force values 
may be different in the positive and negative directions to enable the simulation of nonsymmetric 


envelopes. 


2.3.2 Degrading Model for Potential Weld Failures 


A new hysteresis model was developed incorporating the effects of potential weld fracture 
on the inelastic response of the connection region. Since limited experimental data is available on 


the response of welded connections following a weld failure, a conceptual model was developed 


based on preliminary test data generated at the University of Texas, Austin and the University of 
California, Berkeley. The features of the proposed model are shown qualitatively in Figure 2-5 


and summarized below: 


The primary response is characterized by a bilinear envelope with a yield capacity specified 
by M,. The moment at the instant of weld failure is denoted by M.,._ Presently, this critical 
moment is specified as a function of the yield moment. When additional data becomes available, it 
will be possible to replace the parameter B; by a more comprehensive parameter which reflects 
cumulative damage at the connection. At the onset of weld failure, the primary envelope is 
replaced by a new degraded bilinear representation with reduced stiffness (specified as B2 ki), 
reduced capacity (B; M,) and modified post-yield slope (specified as a function of the initial post- 
yield stiffness, B; kz). Unloading from the new envelope results in a degraded stiffness expressed 
as a function of the new reduced stiffness, B2B4k:. Unloading paths aim the initial stiffness path 
on the negative side, unless the degree of inelasticity causes unloading to reach the post-yield 


stiffness path directly, as demonstated in the loop behavior (Figure 2-5). 


Since weld failures occur primarily in positive bending, it is assumed that the hysteresis 


loops on the negative side will retain the original stiffness and capacity, as shown in Figure 2-5, 


10 


Figure 2-4. Bilinear Hysteretic Model 
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strength 
drop 


(a) Model Parameters 


(b) Loop Behavior 


Figure 2-5. Hysteresis Model for Connection Region 
Before and After Weld Fracture 
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3 MODEL VALIDATION STUDIES 


The models developed in the previous section were implemented in an enhanced version of 
the IDARC computer program. Given the features of this new version, which can model both 
concrete and steel sections simultaneously, a new acronym is used - IDASS, representing Inelastic 
Dynamic Analysis of Structural Systems. Program IDASS was validated using experimental 
results of available component tests and an existing computer program DRAIN-2DX (Prakash et 
al., 1992). The primary purpose of the comparative studies is to validate the new models, the 


details of which are described herein. 


3.1 Member Model Validation 


The new member model developed in Section 2.1, based on distributed flexibility 
formulations, was compared to the two-component model in DRAIN-2DX which uses a 
concentrated plasticity approach. The verification was carried out for two loading conditions: (a) 
- §tatic; and (b) dynamic. In both cases, a sample single-bay single-story structure was used in the 
evaluations, the details of which are shown in Figure 3-1. A post-yield stiffness ratio of 0.05 was 
used in all simulations. Results of the quite analysis are summarized in Table 3-1. The first load 
case corresponds to an elastic response. The second load case produced yielding in the beam 
only. The final load of 45 kips produced yielding in the base of the columns in addition to beam 
yielding. As can be inferred from the Table, the flexibility formulations are identical to those 
based on concentrated plasticity when the state of the element is similar at both ends, viz. elastic 
at both ends or yielding at both ends. Model 2 yields slightly higher values than Model 1. The 
insignificant difference is due to the fact that yielding progresses to only 0.5% across the column 
dimension. Using a predefined hinge length value of about 10% produced results much closer to 
that of the concentrated plasticity model, however, no correlation between the distributed 


flexibility and concentrated plasticity models based on hinge length could be established. 
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120" 


| 180" | 
Section Properties 


Element (EA) (El) My 
Columns 6E05 7.5E06 1,800 
Beam rigid 5.0E06 650 


Figure 3-1. Sample Frame Structure Used in Member Model Validation Studies 
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Table 3-1. IDASS vs. DRAIN-2DX Comparison for Static Loads 


“DRAIN-2DX_ 
IDASS-1 
IDASS-2 


DRAIN-2DX 
IDASS-1 
IDASS-2 

DRAIN-2DX 
IDASS-1 
IDASS-2 


* IDASS - 1 : Model 1, Linear Flexibility Variation 
Model 2, Finite Hinge Length Based on End Moment 


The second phase of numerical testing involved dynamic analysis of the frame subjected to 
seismic loads. A floor weight of 200 kips and zero damping were assumed for the dynamic 
evaluation which resulted in a fundamental period of 0.6 seconds. As in the previous case, three 
load cases were investigated. The 1940 El Centro acceleration record was used as input. The 
seismic evaluations were all carried out at a time step of 0.02 seconds which corresponds to the 
input time step of the acceleration record. The first loading event with a PGA of 0.04g produced 
an elastic response. The next time history at a PGA of 0.08 g produced yielding in the beam only. 
The time history responses obtained using IDASS were identical to those obtained using DRAIN- 
2DX for both cases. Figure 3-2 shows the comparison for the elastic response only. Yielding in 
both columns and beams was observed at a PGA of 0.12g. The comparative response for this 
loading is shown in Figure 3-3. There are slight discrepancies in the response using the two 


programs, however, the difference at maximum amplitude is less than 0.5%. 
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Figure 3-3. Inelastic Response Comparison for Seismic Input (El Centro 0.12g) 
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3.2 Joint Panel Verification 


The joint panel model was verified using two separate sets of available results. The first 
was based on data reported by Charney and Johnson (undated) who used both simplified 
analytical expressions and detailed finite element analysis to compute lateral story displacements 
including the effect of joint panel deformations. In their study, they compared three different 
modeling approaches: (1) element lengths based on center-line dimensions; (2) element lengths 
based on clear span dimensions wherein the joint is considered to be a rigid zone; and (3) element 


lengths based on clear span wherein the joint is modeled as a panel element. 


The results reported in Charney and Johnson’s paper are based on a beam-column 
subassemblage that simulates an interior bay in a moment-resisting frame structure. The basic 
scheme of the beam-column assemblage is shown in Figure 3-4. Equal and opposite loads are 
applied at the beam ends as shown. Results of the analyses using the panel model presented in 
Section 2.2 and implemented in IDASS are tabulated in Table 3-2 and compared to the FEM 
solution reported in Charney and Johnson’s paper. Additionally, a commercial computer program 


_ STAAD-III was used to verify the center-line and rigid joint models. 


This phase of the evaluation was limited to comparing the elastic response of the joint 
panel region. In modeling the joint panel region in IDASS, the cross-section properties assigned 
to the panel were based on the depth of the beam, the depth of the column and the thickness of 


the column web. 


le 


1000 kips 


column span 


1000 kips 


Pp _____——_] 


beam span 


Figure 3-4. Beam-Column Assemblage Used for Panel Model Validation 


18 


Table 3-2. Validation of Joint Panel Model - 


Elastic Interstory Drift Response 


A. Beam W14X132 | FEM Analysis 
Column W21X101 STAAD-III 
IDASS 


B. Beam W14X426| FEM Analysis 

Column W27X178 STAAD-III 
IDASS 

C. Beam W14X426 | FEM Analysis 

Column W36X300 STAAD-III 
IDASS 


drift response values in inches 


* Specimen A: Column height = 240” ; Beam span = 150” 
Specimen B & C: Column height = 360” ; Beam span = 150” 
+ FEM results based on data reported in Chamey and Johnson 


The second phase of verification studies was directed towards computing the inelastic 
response of the joint panel region. To accomplish this, experimental results from a series of tests 
conducted at Lehigh University (Sarkisian, 1985) were used. The general layout of a typical 
specimen is shown in Figure 3-5. W24X62 beam sections and W14X90 column sections were 
used in the testing. The beam-column assemblage was fixed at the base and generally free to 
translate and rotate at the top. Loading was applied by means of hydraulic jacks in equal and 
Opposite directions at each end of the beam. Instrumentation of the panel region provided a direct 


measure of the panel zone deformation. 


19 


3'-6" 


je \ 
W24x 62 


37g"x 34" 


Doubler Piates Connection Plates 


v2'x4" With 7 
Spee. bi. i NO. falelty ASTM-A325 Bolts 
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Full Penetration 
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3  w2' Grade 50 
3/8" Fillet Welds 


4 3/4" Grade 36 
3/8" Fillet Welds 


Figure 3-5. Details of Beam-Column Specimen Used for Validating 
Inelastic Panel Zone Model 


20 


The panel model implemented in program IDASS provides a measure of the change in 
angle (from the original right angle) which can be transformed into panel zone deformation using 


the following relationship: 


(3-1) 


where: 


fede ed. onl a=aldesd,)—2:d.d,.cos0 (3-2) 


d, and d, are the depths of the column and beam, respectively and @ is the change in angle 


between the column and beam rotations at the joint. 


Results of the simulation for two complete cycles are shown in Figure 3-6. Also shown is 
the experimentally observed response. Given the limitations of the bilinear hysteretic model, the 
proposed panel model is capable of representing inelastic shear distortions in the joint with 


acceptable accuracy. 
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(B) ANALYTICAL SIMULATION 


Figure 3-6. Observed vs. Simulated Inelastic Response of Panel Zone 


2g 


3.3. Validation of Hysteresis Model 


The final task in the validation phase of this research effort was the verification of the new 
hysteresis model for the connection region following a sudden weld failure. Preliminary results 
from an experimental study recently conducted at the University of California, Berkeley as part of 


the SAC Joint Venture, were used for this purpose. 


The test setup, specimen details and the computer model used in IDASS for the nonlinear 
evaluation are all shown in Figure 3-7. Three full-scale tests were conducted on a typical beam- 
column connection built according to industry standards prior to the Northridge earthquake. Of 
these, one specimen that was subjected to one full cycle following weld failure was selected for 
the purpose of validating the hysteretis model presented in Section 2-3. Results of the analysis 


are plotted along with the observed experimental response in Figure 3-8. 


It is seen that the proposed model is capable of reproducing experimentally observed 
~ behavior of the welded connection region before and after weld failure. In the present analysis, 
the failure of the weld was specified based on observed experimental data. In an actual analysis of 
a steel frame structure, it must be possible to specify this critical failure point based on separate 


analysis or the use of a cumulative damage model. 
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134 In. 


W36 x 150 (A36) 
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(A) TEST SETUP 
es 134" | 
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(B) IDASS MODEL 


Figure 3-7. Details of Specimen Used For Hysteretic Model Verification 
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Figure 3-8. Comparison of Observed and Simulated Hysteretic Response of 
Connection Region Before and After Weld Fracture 
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4 CONCLUDING REMARKS 


Three primary tasks were undertaken in this project with the aim of developing suitable 
modeling schemes that could be used to analyse the nonlinear dynamic response of SMRFs before 


and after the failure of welded connections in critical regions. 


First, a new member model based on concepts of distributed flexibility was developed. It 
was shown that the new model reproduces results predicted by concentrated plasticity for the case 
of both ends of a member having the same state (elastic or yield conditions). For the case that 
only one end-section of a member yields while the other remains elastic, the flexibility-based 
model produces values that are less than those predicted by the concentrated plasticity model. 
Given the fact that concentrated plasticity always over-predicts observed response, it can be 


concluded that the proposed model may be a better representation of actual inelastic behavior. 


A macromodel representation of panel distortion in the joint region of a moment frame 
was developed. Analytical simulations using the model were compared to results obtained by 
rigorous finite element analysis and to observed experimental behavior. It is established that the 


proposed formulation predicts with acceptable accuracy the inelastic behavior of the panel region. 


Finally, a new hysteresis model was developed to simulate the condition of a sudden weld 
failure. The model was derived conceptually from observed experimental response of such 
connections before and after weld failure. Parameters currently assigned to the model can be 
enhanced in future as more data becomes available. The model was validated using available 


experimental data from a series of tests conducted at the University of California, Berkeley. 
All models described and validated in this report have been incorporated in program 


IDASS. The User Manual for the program and the data sets used to reproduce the results 


presented in this report are included in the Appendices. 
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PROGRAM USER GUIDE 
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A _ 
Inelastic Damage Analysis of Structural Systems - Version 1.0 


USER GUIDE 


INPUT FORMAT 


A free format is used to read all input data. Hence, conventional delimiters (commas, blanks) 
may be used to separate data items. Standard FORTRAN variable format is used to distinguish 
integers and floating point numbers. Input data must, therefore, conform to the specified 
variable type. 


NOTE: Provision is made for a line of text between each set of data items. Refer to the sample 
data files accompanying this Manual. No blank lines are to be input. A zero input will result in 
program default values, where applicable. 


TITLE OF PROBLEM 


e TITLE Alpha-numeric title, upto 80 characters. 
CONTROL DATA (SEE FIGURE A-1) 
e USER_TEXT Reference information: upto 80 characters of text 


e NSO,NFR,NCON,NSTL,NPDEL 


NSO = Number of stories 

NFR = Number of frames 

NCON = Number of different concrete material properties ( = 0 for steel structures) 
NSTL = Number of different steel stress-strain envelopes specified 

NPDEL = 0 ; ignore P-Delta effects; = 1; include P-Delta effects 


NOTES: A structure must be decomposed into a Series of parallel frames. Input is required only 
for non-identical frames, denoted here by the integer variable NFR. The entire group of frames 
can be defined using an L-I-J nodal locater system. This concept is shown graphically in Figure 
A-1 for three different examples. In Figure A-la, the four-story building made up of a total of 
four frames is assumed to have two pairs of identical frames, hence, only two of them need be 
input in IDARC (NFR=2). The cantilever beam/column shown in Figure A-1b is defined as a 
single-story structure with one column line. Likewise, the subassemblage shown in Figure Ic is 
defined as a 2-story structure with three column lines. 
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INTERIOR FRAMES (c) Typical Beam-Column 
Subassemblage 
(a) Typical Building 
Figure A-1. Frame Discretization and Nodal Identification 
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ELEMENT TYPES (SEE FIGURE A-1) 
‘e USER_TEXT Reference information: upto 80 characters of text 


e MCOL,MBEM,MWAL,MEDG, MTRN,MSPR,MJNT 


MCOL = No. of types of columns 

MBEM = No. of types of beams 

MWAL = No. of types of shear walls 
MEDG = No. of types of edge columns 
MTRN = No. of types of transverse beams 
MSPR = No. of types of rotational springs 
MJNT = No. of types of joints 


NOTES: Elements are grouped into identical sets based on cross-section data and initial 
conditions such as axial loads. For example, in the interior frame shown in Figure A-la, 
assuming identical interior and exterior columns in each floor, only 8 column types are needed 
to define all 16 elements, i.e., 2 types per each level as shown in the Figure. 


ELEMENT DATA 
e USER_TEXT Reference information: upto 80 characters of text 


e NCOL,NBEM,NWAL,NEDG,NTRN,NSPR,NMR, NJNT 


NCOL = Total number of columns 

NBEM = Total number of beams 

NWAL = Total number of shear walls 
NEDG = Total number of edge columns 
NTRN = Total number of transverse beams 
NSPR = Total number of rotational springs 
NJNT = Total number of joints 

NMR = Total number of moment releases 


NOTES: NMR is used to specify moment releases (hinge locations) at member ends. Releasing 
a moment at a member end results in a hinge condition at that end thereby disallowing moments 
to develop at the section. Moment releases may not be specified at both ends of a member. 


UNIT SYSTEM 
e USER_TEXT Reference information: upto 80 characters of text 
elU System of units 

= 1, inch,kips 

= 2, mm, kN 
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FLOOR ELEVATIONS (SEE FIGURE A-2) 


e USER_TEXT Reference information: upto 80 characters of text 


e HIGT(1),HIGT(2)...HIGT(NSO) Elevation of each story from the base, 
beginning with the first floor level. 


DESCRIPTION OF IDENTICAL FRAMES 


e USER_TEXT Reference information: upto 80 characters of text 
e NDUP(1),NDUP(2)...NDUP(NFR) Number of duplicate frames for each 
of the NFR frames 


NOTES: In the sample structure shown in Figure A-1, there are four frames. However, the two 
interior frames are identical as are the exterior frames. In this case, NFR=2, and NDUP(1) = 
NDUP(2) = 2. 


PLAN FIGURATION 


e USER_TEXT Reference information: upto 80 characters of text 


e NVLN(1),NVLN(2)...NVLN(NFR) Number of column lines (or J-locater points) 
in each frame. 


NOTES: A set of NVLN points for each frame should define completely the column lines 
necessary to specify every vertical element in that frame. If a beam element is subdivided into 
two or more segments, then the number of column lines specified must include these internal 
beam nodes as well. 


NODAL WEIGHTS (SEE FIGURE A-2) 
e USER_TEXT Reference information: upto 80 characters of text 


e LEVEL, IFR(1), WVT(1) . WVT(NVLN(1)) 
IFR(2), WVT(1)...WVT(NVLN(2)) 
hs repeat for NFR frames 


e (repeat upto NSO levels) 


LEVEL = Story level number 
IFR(J) = Frame number 
WVT(K) = Nodal weight 
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HIGT(1) 


FRAME #1 


(numbers shown at nodes = nodal weights) 


INPUT DATA: 3.0, 6.0, 6.0, 6.0, 3.0 
5.0, 0.0, 10.0, 0.0, 5.0 


5.0, 0.0, 10.0, 0.0, 5.0 


Figure A-2. Floor Heights and Nodal Weights 
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HIGT(3) 


ENVELOPE GENERATION OPTION 


e USER_TEXT Reference information: upto 80 characters of text 
e IUSER Code for specification of user properties 
= 0, requires IDARC generated envelopes for 
atleast one element 
= 1, complete moment-curvature envelope data 
to be provided by user 


SKIP THIS INPUT IF IUSER = 1 


e USER_TEXT Reference information: upto 80 characters of text 
e IM,FC,EC,EPSO,FT,EPSU,ZF 


e (repeat for each of the NCON concrete types) 


IM = Concrete type number 

FC = Unconfined compressive strength 

EC = Initial Young's Modulus of concrete 
EPSO = Strain at max. strength of concrete (%) 
FT = Stress at tension cracking 

EPSU = Ultimate strain in compression (%) 

ZF = Parameter defining slope of falling branch 


Default Values: EC =57Vf' ksi; EPSO=0.2% ; FT =0.12*FC; 
EPSU and ZF are computed from cross-section data. 


SKIP THIS INPUT IF IUSER = 1 


e USER_TEXT Reference information: upto 80 characters of text 
e IM,FS,FSU,ES,ESH,EPSH Characteristics of steel stress-strain curve: 


e (repeat for each of the NSTL steel types) 


IM = Steel type number 

FS = Yield strength 

FSU = Ultimate strength 

ES = Modulus of elasticity 

ESH = Modulus of strain hardening 
EPSH = Strain at start of hardening (%) 


Default Values: 
FSU=14*FS; ES =29,000 ksi; ESH = (ES / 60) ksi; EPSH = 3.0% 
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Stress, © 


Ut 


Strain, € 


FT 


Figure A-3. Stress-Strain Curve for Unconfined Concrete 


Strain, € 


Symmetric 


Figure A-4. Stress-Strain Specification for Steel Sections and Reinforcing Bars 
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e USER_TEXT Reference information: upto 80 characters 
e NHYS Number of types of hysteretic rules 


HYSTERETIC MODEL PARAMETERS 


e IR, IHYSTYP, PARAM1, PARAM2, PARAM3, PARAM4, PARAMS, PARAM6 
e (NHYS lines of data) 


( SEE FIGURES A-5, A-6 AND A-7 ) 


IR = Parameter Set Number 

IHYSTYP = 1 , Reinforced Concrete 
2 , Steel (non-degrading) 
3 , Steel (degrading) 


Input values of hysteretic parameters based on the choice of IHYSTYP as follows: 


IHYSTYP = 1 IHYSTYP =2 IHYSTYP = 3 
PARAM 1 Stiffness degrading Post-yield stiffness Post- Yield Stiffness 
coefficient ratio Ratio, @ 


PARAM2 Energy-based strength | Not used, input 0.0 Strength Reduction 
epiveinien decay parameter ohana Ratio after Weld 
Failure, 
strength decay Ratio, B, 
parameter 
closing parameter Reduction Ration, 


Stiffness Ratio, B, 
PARAM6 Not used, input 0.0 Not used, input 0.0 Critical Force Factor 
Failure, B 
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1. Modeling of Stiffness 
Degradation 


c AF 2- Modeling of Strength 
Deterioration 


Feder Fe (AeA Et- A2*D) 


FY Uult 
E = Area of F-u loops 
D = Ductility 
A1 = PARAM2 
A2 = PARAMS 


3. Modeling of Slip 
or Pinching Behavior 


PYP 


PCP 


~ PARAM4*PYP 


Figure A-5. Hysteretic Parameters for RC Sections 
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PARAM 1*k 4 


PARAM1 “Ky EE 


Figure A-6. Bilinear Nondegrading Model 


Figure A-7. Hysteresis Model for Steel Connections Including Weld Fracture 
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e USER_TEXT Reference information: upto 80 characters of text 
e IUCOL Type and option for column section input 


=1; Reinforced Concrete; cross-section data input 

=2; Reinforced Concrete; moment-curvature input 

=3; Steel; cross-section data - bare steel (symmetric) 

=4; Steel; moment-curvature input 

=5; Composite (steel and concrete) and nonsymmetric section 


IF IUCOL = 1, CONTINUE WITH SET E1 

IF IUCOL = 2, GO TO SET E2 

IF IUCOL = 3, GOTO SET E3 

IF IUCOL = 4, GOTO SET E4 (IUCOL=5, unavailable) 


DATA SET El 


e USER_TEXT Reference information: upto 80 characters of text 


For each column type, input the following: 


e ICTYPE Type of column 
= 1; rectangular (DEFAULT) 
= 2; circular 
Rectangul ion Data: (Figure A-8) 
General data: 


e KC,IMC,IMS,AN,AMLC,RAMC1,RAMC2 
Bottom section: 
e KHYSC, D, B, DC, AT, HBD, HBS, CEF 
Top section (skip if symmetric, see note below): 
e KHYSC, D, B, DC, AT, HBD, HBS, CEF 


NOTE: If KHYSC for bottom section is input with negative sign, section is assumed to have 
identical properties for bottom and top section; no input is required for top section 


KC = Column type number 

IMC = Concrete type number 

IMS = Steel type number 

AN = Axial load 

AMLC = Center-to-center column height 
RAMC1 = Rigid zone length at bottom 
RAMC2 = Rigid zone length at top 


AMLC 


Note moment 
sign convention 


Typical Column Cross-Section 


Typical Column Line 


O} @ 


Minimal Confinement Nominal Confinement Well Confined 
CEFF = 0.5 CEFF = 0.66 CEFF = 1.0 


Effectiveness of Confinement for Some Typical 
Hoop Arrangements 


Figure A-8. Rectangular Concrete Column Input Details 
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KHYSC = Hysteretic rule number (may be negative) 
D = _ Depth of column 
B = Width of column 
= _ Distance from centroid of reinforcement to face of column 
AT = Area of reinforcement on one face 
HBD =Hoop bar diameter 
HBS = Hoop bar spacing 
CEF = Effectiveness of column confinement 


Return to input of ICTYPE. When done, go to SET F 


Circular Section Data: (Figure A-9) 
General data: 
e KC,IMC,IMS,AN,AMLC,RAMC1,RAMC2 
e KHYSC, AN,DO,CVR,DST,NBAR,BDIA,HBD,HBS 


KC =Colum Type number 

IMC = Concrete type number 

IMS = Steel type number 

AMLC = Center-to-center column height 
RAMC1 = Rigid arm bottom 

RAMC2 = Rigid arm top 


KHYSC = Hysteretic Rule number 

AN = Axial load on the column 

DO = Outer diameter of column 

CVR = Cover to center of hoop bar 

DST = Distance between centers of long. bars 
NBAR = Number of longitudinal bars 

BDIA = Diameter of longitudinal bar 

HBD = Diameter of hoop bar 

HBS = Spacing of hoop bars 


Return to input of ICTYPE. When done go to SET F. 
SET E2: REINFORCED CONCRETE - MOMENT CURVATURE INPUT (Figure A-10) 


e USER_TEXT Reference information: upto 80 characters of text 
General Data: 

e KC, AMLC, RAMC1, RAMC2 

Bottom section: 

e KHYSC, EI,EA,GA, PCP,PYP,U YP,UUP,EI3P,PCN,PYN,UYN,UUN,EI3N 


Top section (skip if symmetric, see note below): 
e KHYSC, EI,EA,GA, PCP,PYP,UYP,UUP,EI3P,PCN,PYN,UYN,UUN,EBN 


Note: A negative sign for KHYSC for bottom Section indicates similar properties for top section. 
e (repeat for each of MCOL sections) 
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a 


CVR | 
DO 


CROSS-SECTION 


CIRCULAR HOOPS SPIRAL HOOPS 


Figure A-9. Circular Concrete Column Input Details 


Force EISPor GASP 


Deformation 
UUP 


Note: Force = Moment or Shear 
Deformation = Curvature, Rotatior 


or Strain 
EI3Nor 


Figure A-10. Moment-Curvature Input for RC Sections 


A-14 


KC = Column type number 

-AMLC = Column Length 

RAMCI = Rigid Arm (Bottom) 

RAMC2 = Rigid Arm (Top) 

KHYSC = Hysteretic rule number (may be negative) 
EI = Initial Flexural Rigidity (EI) 

EA = Axial Stiffness (EA/L) 

GA = Shear Stiffness (Shear modulus*Shear Area) 
PCP = Cracking Moment (positive) 

PYP = Yield Moment (positive) 

UYP = Yield Curvature (positive) 

UUP = Ultimate Curvature (positive) 

EI3P = Post yield Flexural Stiffness (positive) 
PCN = Cracking Moment (negative) 

PYN = Yield Moment (negative) 

UYN = Yield Curvature (negative) 

UUN = Ultimate Curvature (negative) 

EI3N = Post yield Flexural Stiffness (negative) 


ET E3: STEEL - CR -SECTION INPUT (FIGURE A-11) 


e USER_TEXT Reference information: upto 80 characters of text 


Section Data: 
e KC, IMS, AMLC, RAMC1, RAMC2, AN, D, BF, TF, TW, AX, AY, IZ, SX, ZX 
e (repeat for each of MCOL sections) 


KC = Column type number 

IMS = Steel stress-strain property number 
AMLC = Column Length 

RAMC1 = Rigid zone length (Bottom) 
RAMC2 = Rigid zone length (Top) 
AN = Axial load 

D = Total depth of section 

BF = Flange width 

TF = Flange thickness 

TW = Web thickness 

AX = Cross-sectional area 

AY = Shear Area 

IZ = Moment of Inertia 

SX = Elastic Section Modulus 

ZX = Plastic Section Modulus 


NOTE: Zero inputs for D, BF, TF or TW require non-zero inputs for AX, IZ, SX and ZX. 
Zero inputs for AX, IZ, SX or ZX require non-zero inputs for D, BF, TF and TW. 
Shear deformations will be ignored if AY = 0 
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4 UE Be ee 


rf | 


Figure A-11. Input Parameters for Symmetric Steel W-Sections 


AT 1 


Figure A-12. Input Details for RC Beam Section 
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‘SET E4; STEEL - MOMENT CURVATURE INPUT 


e USER_TEXT Reference information: upto 80 characters of text 


General Data: 

e KC, AMLC, RAMC1, RAMC2 

Bottom section: 

e KHYSC, EI, EA, GA, PYP, PYN 

Top section (skip if symmetric, see note below): 
e KHYSC, EI, EA, GA, PYP, PYN 


Note: If KHYSC for bottom section is input with negative sign, section is assumed to have 
identical properties for top section; skip top section input 


e (repeat for each of MCOL Sections) 


KC = Column type number 

AMLC = Column Length 

RAMC1 = Rigid zone (Bottom) 

RAMC2 = Rigid zone (Top) 

KHYSC = Hysteretic rule number (may be negative) 
EI = Initial flexural Rigidity (EI) 

EA = Axial stiffness (EA/L) 

GA = Shear stiffness (Shear modulus*Shear Area) 
PYP = Yield moment (positive) 

PYN = Yield moment (negative) 


SKIP THIS INPUT IF THE STRUCTURE HAS NO BEAMS (NBEM=0) 


e USER_TEXT Reference information: upto 80 characters of text 
e IUBEM Type and option for beam section input 


=1; Reinforced Concrete; cross-section data input 

=2; Reinforced Concrete; moment-curvature input 

= 3; Steel; cross-section data (bare steel, symmetric) 

= 4; Steel; moment-curvature input 

=5; Composite (steel and concrete) and nonsymmetric section 


IF TIUBEM = 1, CONTINUE WITH SET F1 

IF IUBEM = 2, GOTO SET F2 

IF IUBEM = 3, GOTO SET F3 

IF IUBEM = 4, GOTO SET F4 (IUBEM =5, unavailable) 
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DATA SET F1 (Figure A-12) 


e USER_TEXT Reference information: upto 80 characters of text 


General data: 
e KB,IMC,IMS,AMLB,RAMB1,RAMB2, 
Left section: 


e KHYSB, D, B, BSL TSL, BC, AT1, AT2, HBD, HBS 
Right section (skip, if symmetric, see note below): 
e KHYSB, D, B, BSL TSL, BC, AT1, AT2, HBD, HBS 


e (repeat for each of MBEM sections) 


Note: If KHYSB for left section is input with negative sign, section is symmetric and input for 
right section should be omitted. 


KB = Beam type number 

IMC = Concrete type number 

IMS = Steel type number 

AMLB = Member length 

RAMB1 = Rigid zone length (left) 

RAMB2 = Rigid zone length (right) 

KHYSB = Hysteretic rule number (may be negative) 
D = Overall depth 

B =Lower width 

BSL = Effective slab width (=B for rectangular section) 
TSL = Slab thickness (= 0 for rectangular section) 
BC = Cover to centroid of steel 

AT1 = Area of bottom bars 

AT2 = Area of top bars 

HBD = Diameter of stirrup bars 

HBS = Spacing of stirrups 


SET F2: REINFORCED CONCRETE - MOMENT CURVATURE INPUT (Figure A-10) 


e USER_TEXT Reference information: upto 80 characters of text 


General Data: 

e KB, AMLB, RAMB1, RAMB2 

Left section: 

e KHYSB, EI,GA, PCP,PYP,UYP,UUP,EI3P,PCN,PYN,UYN,UUN,EI3N 
Right section (skip if symmetric, see note below): 

e KHYSB, EI,EA,GA, PCP,PYP,UYP,UUP,EI3P,PCN,PYN,UYN,UUN,EI3N 


Note: If KHYSB for left section is input with negative sign, section is assumed to be symmetric, 
and right section data input should be omitted.. 


e (repeat for each of MBEM sections) 


_KB = Beam type number 

AMLB = Beam length 

RAMBI1 = Rigid zone (left) 

RAMB2 = Rigid zone (right) 

KHYSB = Hysteretic rule number (may be negative) 
EI = Initial Flexural Rigidity (EI) 

GA = Shear Stiffness (Shear modulus*Shear Area) 
PCP = Cracking Moment (positive) 

PYP = Yield Moment (positive) 

UYP = Yield Curvature (positive) 

UUP = Ultimate Curvature (positive) 

EI3P = Post yield Flexural Stiffness (positive) 
PCN = Cracking Moment (negative) 

PYN = Yield Moment (negative) 

UYN = Yield Curvature (negative) 

UUN = Ultimate Curvature (negative) 

EI3N = Post yield Flexural Stiffness (negative) 


ET F3: STEEL - CR -SECTION INPUT 


e USER_TEXT Reference information: upto 80 characters of text 


Section Data: 


e KB, IMS, AMLC, RAMC1, RAMC2, D, BF, TF, TW, AX, AY, IZ, SX, ZX 
e (repeat for each of MBEM sections) 


KC = Column type number 

IMS = Steel stress-strain property number 
AMLC = Column Length 

RAMCI1 = Rigid zone length (Bottom) 
RAMC2 = Rigid zone length (Top) 
AN = Axial load 

D = Total depth of section 

BF = Flange width 

TF = Flange thickness 

TW = Web thickness 

AX = Cross-sectional area 

AY = Shear Area 

IZ = Moment of Inertia 

SX = Elastic Section Modulus 

ZX = Plastic Section Modulus 


SEE NOTES FOR SET E-3 
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SET F4: STEEL _- MOMENT CURVATURE INPUT 


e USER_TEXT Reference information: upto 80 characters of text 


General Data: 

e KB, AMLB, RAMB1, RAMB2 

Left section: 

e KHYSB, EI, GA, PYP, PYN 

Right section (skip if symmetric, see note below): 
e KHYSB, EI, GA, PYP, PYN 


Note: If KHYSB for left section is input with negative sign, section is assumed to be symmetric, 
and right section data input should be omitted.. 


e (repeat for each of MBEM sections) 


KB = Beam type number 

AMLB = Beam Length 

RAMBI = Rigid zone (left) 

RAMB2 = Rigid zone (right) 

KHYSB = Hysteretic rule number (may be negative) 
EI = Initial flexural Rigidity (ED) 

GA = Shear stiffness (Shear modulus*Shear Area) 
PYP = Yield moment (positive) 

PYN = Yield moment (negative) 


e USER_TEXT Reference information: upto 80 characters of text 


e IUWAL Type of wall input 
= (; Cross-section input 
= 1; Moment-curvature and shear-strain input 


IF IUWAL = 1, GO TO SET G2 


ET Gl: CR -SECTION INPUT 
e USER_TEXT Reference information: upto 80 characters of text 
General Data: 


e KW,IMC,KHYSW(1),KHYSW(2),KHYSW(3),AN,AMLW,NSECT 
For each of the NSECT Sections, input the following 

e KS, IMS,DWAL,BWAL,PT,PW 

e repeat NSECT times 
e repeat for each of MWAL sections 
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SECTION 1 SECTION 2 SECTION 3 


. eters by SAE AE) Cae 
BWAL(1) BWAL(3) 
ene! wine 


DWAL(1) DWAL(2) DWAL(3) (section with 
edge columns) 


} 
Ae (section without 


edge columns) 


BWAL(1)=BWAL(2)=BWAL(3) 


DWAL(1) DWAL(2) DWAL(3) 


- ia 


pone 
JSS ee ela 


:— SECTION AA 


Figure A-13. Concrete Shear Wall Input Details 
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KS = Section number 

‘IMS = Steel type number 

DWAL = Depth of section 

BWAL = Width of section 

PT = Vertical reinforcement ratio (%) 
PW = Horizontal reinf ratio (%) 


KW _ = Shear wall type number 

IMC = Concrete type number 

KHYSW(1) = Hysteretic Rule Number (bottom) 
KHYSW(2) = Hysteretic Rule Number (top) 
KHYSW(3) = Hysteretic Rule Number (shear) 
AN = Axial load 

AMLW = Height of shear wall 

NSECT = Number of Sections 


SET G2; MOMENT CURVATURE INPUT (Figure A-10) 


e USER_TEXT Reference information: upto 80 characters of text 


General Data: 

e KW, AMLW, EAW 

Flexure - Bottom section: 

e KHYSW, EI,PCP,PYP,UYP,UUP,EI3P, PCN,PYN,UYN,UUN,EI3N 
Flexure - Top section (skip if symmetric, see note below) 

e KHYSW, EI,PCP,PYP,UYP,UUP,EI3P, PCN,PYN,UYN,UUN,EI3N 
Shear properties: 

e KHYSW, GA,PCP,PYP,UYP,UUP,GA3P, PCN,PYN,UYN,UUN,GA3N 


Note: If KHYSW for bottom section is input with negative sign, section is symmetric, hence, do 
not input top section data 
e repeat for each of MWAL sections 


Flexural data: 

KW = Wall type number 

AMLW = Wall length 

EAW = Axial Stiffness (EA/L) 

KHYSW = Hysteretic rule number (may be negative) 
EI = Initial flexural stiffness (ED) 

PCP = Cracking Moment (positive) 

PYP = Yield Moment (positive) 

UYP = Yield Curvature (positive) 

UUP = Ultimate Curvature (positive) 

EI3P = Post Yield Flexural Stiffness (positive) 
PCN = Cracking Moment (negative) 
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PYN = Yield Moment (negative) 
UYN = Yield Curvature (negative) 
UUN = Ultimate Curvature (negative) 
_EI3N = Post yield Flexural Stiffness (negative) 


Shear data: 

KHYSW = Hysteretic Rule Number 

GA = Initial Shear Stiffness (shear modulus*area) 
PCP = Cracking Shear (positive) 

PYP = Yield Shear (positive) 

UYP = Yield Shear strain (positive) 

UUP = Ultimate Shear strain (positive) 
GA3P =Post Yield Shear Stiffness (positive) 
PCN = Cracking Shear (negative) 

PYN = Yield Shear (negative) 

UYN = Yield Shear strain (negative) 

UUN = Ultimate Shear strain (negative) 
GA3N =Post Yield Shear Stiffness (negative) 


SKIP THIS INPUT IF THE STRUCTURE HAS NO EDGE COLUMNS 
Do not duplicate edge column data if already input as part of shear wall section 


e USER_TEXT Reference information: upto 80 characters of text 
e KE,IMC,IMS,AN,DC,BC,AG,AMLE,ARME 


KE = Edge column type number 
IMC = Concrete type number 
IMS = Steel type number 

AN = Axial load 

DC = Depth of edge column 

BC = Width of edge column 
AG =Gross area of main bars 
AMLE = Member length 

ARME = Arm length 


THIS INPUT NOT REQUIRED IF STRUCTURE HAS NO TRANSVERSE BEAMS 
e USER_TEXT Reference information: upto 80 characters of text 


eKT,AKV,ARV,ALV 
e (repeat for each of MTRN types) 
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FRAME 1 
(1=1) 


WALL 


S aeeeieee 
[t 
D 
TRANSVERSE 
BEAM Waraee 
b 


WALL 3 3 
ARV = 0.33(tB+bD)G 


Column Line 1 Column Line 2 
(J=1) (J=2) 


HORIZONTAL GROUND MOTION 


Figure A-14. Transverse Beam Input Parameters 


KT = Transverse beam type number 
AKV'= Vertical Stiffness 

ARV = Torsional Stiffness 

ALV = Arm length 


NOTES: 


1. Transverse elements are assumed to remain elastic. The degree of fixity at the ends 
will depend on the state(cracked/yielded) of the joint and the members that frame into 
the joint before and during the application of load. If the entire region is expected to 
stay elastic, then the vertical stiffness should be computed as AKV = 12 El/L’. In the 
extreme case that one of ends do not transmit stiffness due to yielding of adjoining 
members or deterioration of the joint, then AKV =3 EI/L’.. An intermediate value is 
a good average approximation. 

2. If duplicate frames are present, extreme care should be taken in specifying 
transverse beam properties. The program multiplies the input values by the number of 
duplicate frames to which they are attached. For example, for the frames shown in 
Figure A-1, NDUP(1) = NDUP(2) = 2. The program will factor the input stiffness 
values by (NDUP(1)*NDUP(2))=4.0. Input stiffnesses should, therefore, be modified 
to account for this effect. If the modeling of transverse elements is crucial to the 
analysis, the use of duplicate frames should be avoided. 
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THIS INPUT NOT REQUIRED IF ROTATIONAL SPRINGS ARE NOT SPECIFIED 


e KHYSR, EI,PCP,PYP,UYP,UUP,EI3P, PCN,PYN,UYN,UUN,EI3N 
e (repeat for each of MSPR springs) 


KHYSR = Hysteretic Rule Number 

EI = Initial Rotational Stiffness 

PCP = Cracking moment (positive) 

PYP = Yield moment (positive) 

UYP = Yield rotation (positive, radians) 
UUP = Ultimate rotation (positive, radians) 
EI3P = Post-yield stiffness ratio (positive) 
PCN = Cracking moment (negative) 

PYN = Yield moment (negative) 

UYN = Yield rotation (negative) 

UUN = Ultimate rotation capacity (negative) 
EI3N = Post yield stiffness ratio (negative) 


NOTES: Spring properties, unlike other element types, are specified in terms of moment and 
rotation (in radians). The envelope follows the same nonsymmetric trilinear pattern as 
Shown in Figure A-10. 


e USER_TEXT Reference information: upto 80 characters of text 
e KJ, KHYSJ, DP, BP, TP, G, PYP, PYN > 


KJ = Panel type number 

KHYSJ = Hysteretic rule number (may be negative) 

DP = Depth of panel (typically the depth of the beam) 

BP = Width of panel (typically the depth of the column) 

TP = Thickness of Panel (typically the thickness of the column web) 
G = Shear modulus 

PYP = Yield shear (positive) 

PYN = Yield shear (negative) 
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NOTE: Element connectivity is established through the 3 positional locaters described in Figure 


A-1: a story level, a frame number and a column line. The hypothetical structure shown 
below is used to demonstrate the input format. Only a representative data set is shown. 
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Figure A-15. Element Connectivity 
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SKIP THIS INPUT IF THE STRUCTURE HAS NO COLUMNS 


e USER_TEXT Reference information: upto 80 characters of text 
e MITC ICJC,LBC,LTC 
e (NCOL lines of data) 


M =Column number 

ITC = Column type number 

IC = Frame number 

JC = Column Line number 

LBC = Story level at bottom of column 
LTC = Story level at top of column 


SKIP THIS INPUT IF STRUCTURE HAS NO BEAMS 


e USER_TEXT Reference information: upto 80 characters of text 


e M,ITB,LB,IB,JLB,JRB 
e (NBEM lines of data) 


M = Beam number 

ITB = Beam type number 

LB = Story level 

IB = Frame number 

JLB = Column Line number of left section 
JRB = Column Line number of right section 


e USER_TEXT Reference information: upto 80 characters of text 
e M,ITW,IW,JW,LBW,LTW 
e (NWAL lines of data) 


M = Wall number 

ITW = Wall type number 

IW = Frame number 

JW =Column line number 
LBW = Story level at bottom 
LTW = Story level at top 


e USER_TEXT Reference information: upto 80 characters of text 
e MITE, IE,JE,LBE,LTE 
e (NEDG lines of data) 


M = Edge column number 

ITE = Edge column type number 

IE = Frame number 

JE = Column line number 

LBE = Story level at bottom of column 
LTE = Story level at top of column 


SKIP THIS INPUT IF STRUCTURE HAS NO TRANSVERSE BEAMS 


e USER_TEXT Reference information: upto 80 characters of text 
e M,ITT,LT,.IWT,JWT,IFT,JFT 
e (NTRN lines of data) 


M = Transverse beam number 

ITT = Transverse beam type number 

LT = Story level 

IWT = Frame number of origin of transverse beam* 
JWT = Column line of origin of transverse beam* 
IFT = Frame number of connecting wall or column 
JFT = Column line of connecting wall or column 


NOTES: *For beam-to-wall connections, IWT and JWT refer to the I,J locations of the wall. 


SKIP THIS INPUT IF ROTATIONAL SPRINGS ARE NOT SPECIFIED 


e USER_TEXT Reference information: upto 80 characters of text 
e M, ISP, JSP, LSP, KSPL 
e (NSPR lines of data) 


M = Spring number 

ISP = Frame number 

JSP = Column line number 

LSP = Story level 

KSPL = Relative spring location as follows: 
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Code for KSPL-> =1, spring on beam, left of joint 
= 2, spring on column, top of joint 
= 3, spring on beam, right of joint 
= 4, spring on column, bottom of joint 


NOTE: The number of springs at a joint is limited to one less than the total number of members 
framing into the joint 


KSPL = 1 ee 


KSPL =2 KSPL =4 


Figure A-16. Spring Location Specification 


© 


e USER_TEXT Reference information: upto 80 characters of text 
e IJ, ITJ, IFJ, JIT, LIT 
e (NJNT lines of data) 


IJ = Joint panel number 
ITJ = Joint panel type 

IFJ = Frame number 

JJT = Column Line number 
LJT = Story level 
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SKIP THIS INPUT IF MOMENT RELEASES ARE NOT REQUIRED (NMR = 0) 


e USER_TEXT Reference information: upto 80 characters of text 
e IDM, IHTY, INUM, IREG 
e (NMR lines of data) 
IDM = ID number 
IHTY = Element type using following code 
CODE: 1 = COLUMN 

2 =BEAM 

3 = WALL 
INUM = Column, Beam or Wall number 
IREG = Location of hinge or moment release 


= 1, BOTTOM or LEFT 
= 2, TOP or RIGHT 


e USER_TEXT Reference information: upto 80 characters of text 
e IOPT Option for continuing analysis 
= 0 , STOP (Data check mode) 
= 1, Inelastic incremental analysis with static loads 
= 2, Monotonic "pushover" analysis including static loads (if specified) 
= 3 , Inelastic dynamic analysis including static loads (if specified) 
= 4 , Quasi-static cyclic analysis including static loads (if specified) 


Notes: It is generally advisable to use the "data check" mode for the first trial run of a new 
data set. The program performs only minimal checking of input data. Structural 
elevation plots generated by IDARC help identify errors in connectivity specification. 
Since IDARC prints all input data almost immediately after they are read, the task of 
detecting the source of input errors is generally expedited. It is also important to verify 
all printed output, before carrying out a time-history analysis. 


OPTION I permits an independent nonlinear static analysis. Static loads are input in 
data set T1. OPTIONS 2 - 4 may be combined with long-term static loads which is 
input in data set T1. Initial forces and moments generated by the static loads will 
remain on the structure for all the other options. If a static analysis is not performed, 
the axial loads input as part of column properties will be used as initial axial forces. 


NOTE: THIS INPUT IS REQUIRED FOR ALL ANALYSIS OPTIONS. 
Control Information 


e USER_TEXT Reference information: upto 80 characters of text 
e NLU,NLJ,NLM,NLC 


NLU = No. of uniformly loaded beams 
NLJ = No. of laterally loaded joint 

NLM = No. of specified nodal moments 
NLC = No. of concentrated vertical loads 


IF NLU = NLJ = NLM = NLC = 0, and IOPT = 2, CONTINUE TO SET T2. 


IF NLU = LJ = NLM = NLC =0, and IOPT = 3, CONTINUE TO SET T3. 
IF NLU = \\LJ =NLM = NLC = 0, and IOPT = 4, CONTINUE TO SET T4. 


Next Data Set: 
e JSTP,IOCRL 


JSTP = No. of incremental steps in which to apply the static loads (default = 1 step) 
IOCRL = Steps between printing output (If IOCRL=0, only final results will be printed) 


NOTES: Dead and live loads that exist prior to the application of seismic or quasi-static cyclic 
loads can be input in this section. Such loads are typically specified through uniformly 
loaded beam members. An option is also available for lateral load analysis and the 
specification of nodal loads at joints. When used in conjenction with Options 2-4, the 
resulting forces are carried forward to the monotonic, dynamic and quasi-static 
analysis. 


Uniformly Loaded Beam Data 


SKIP THIS INPUT SECTION IF NLU=0 


e USER_TEXT Reference information: upto 80 characters of text 
e IL, IBN, FU 
e (NLU lines of data) 


IL = Load number 
IBN = Beam number 
FU = Magnitude of load (Force/length) 


Laterally Loaded Joints 
SKIP THIS INPUT SECTION IF NLJ=0 


e USER_TEXT Reference information: upto 80 characters of text 
e IL, LF, IF,,FL 
e (NLJ lines of data) 


IL = Load number 

LF = Story level number 
IF = Frame number 

FL = Magnitude of load 


Nodal Moment Data 
SKIP THIS INPUT SECTION IF NLM=0 


e USER_TEXT Reference information: upto 80 characters of text 
e IL, IBM, FM1, FM2 
e (NLM lines of data) 


IL = Load number 

IBM = Beam number 

FM1 = Nodal moment (left) (See Figure A-9 for beam moment sign convention) 
FM2 = Nodal moment (right) 


Concentrated Vertical Loads 
SKIP THIS INPUT SECTION IF NLC=0 


e USER_TEXT Reference information: upto 80 characters of text 
e IL, IFV, LV, JV, FV 
e (NLC lines of data) 


IL = Load number 

IFV = Frame number 

LV = Story level number 
JV = Column line number 
FV = Magnitude of load 


IF IOPT = 2, CONTINUE TO SET T2. 
IF IOPT = 3, CONTINUE TO SET T3. 
IF IOPT = 4, CONTINUE TO SET T4. 
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e USER_TEXT Reference information: upto 80 characters of text 
e PMAX, MSTEPS 


PMAX = Estimate of base shear strength coefficient (ratio of lateral load capacity to total weight) 
MSTEPS = Number of steps in which to apply the monotonically increasing load 


DEFAULT VALUES: PMAX = I1/NSO + 0.01*NSO ; MSTEPS = 40 


NOTES: The program uses the PMAX value only to determine the load steps for the push-over 
analysis. The prescribed base shear (product of PMAX and total structure weight) is 
applied incrementally in MSTEPS steps as an inverted triangular load, until the top 
story displacement reaches 2% of the total structure height OR the specified PMAX is 
reached. If the program output shows a linear shear vs. deformation plot, the base 
Shear estimate is too low. If the maximum displacement is reached too quickly 
(indicated by too few points in the plot), the estimate is too high. 


IF IOPT = 2 , STOP HERE 


e USER_TEXT Reference information: upto 80 characters of text 
e GMAXH,GMAXV,DTCAL,TDUR,DAMP 


GMAXH = Peak horizontal acceleration (g's) 
GMAXV = Peak vertical acceleration (g's) 
DTCAL = Time step for response analysis (secs) 
TDUR = Total duration of analysis (secs) 
DAMP = Damping coefficient (% of critical) 


NOTES: The input accelerogram is scaled uniformly to achieve the specified peak acceleration. 
DTCAL should not exceed the time interval of the input wave, DTINP. The ratio 
(DTINP/DTCAL) must yield an integer number. TDUR may be less than the total 
duration of the earthquake. If TDUR is greater than the total time duration of the 
input wave, a free vibration analysis of the system will result for the remaining time. 


INPUT WAVE DATA 


e USER_TEXT Reference information: upto 80 characters of text 
e IWV,NDATA,DTINP 
IWV =O, Vertical component of acceleration not included 

= 1, Vertical component of acceleration is included 


NDATA = Number of points in earthquake wave files 
DTINP = Time interval of input wave 
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WAVE TITLE 


e NAMEW Alpha-numeric title for input wave upto 80 characters 


FILENAME - HORIZONTAL COMPONENT 
e WHFILE 
Name of file (with extension) from which to read horizontal component of earthquake record 


Note: Filename should not exceed 12 characters incl. extension 


e WINPH(D),I=1,NDATA ; 
Horizontal component of earthquake wave (NDATA points) is read from the file WHF ILE 


WAVE DATA - VERTICAL COMPONENT 

SKIP THIS INPUT IF IWV .EQ.0 

e WVFILE 

Name of file (with extension) from which to read vertical component of earthquake record 
Note: Filename should not exceed 12 characters 


e WINPV(I),I=1,NDATA 
Vertical component of earthquake wave (NDATA points) read from the file WVFILE 


NOTES: Accelerogram data may be input in any system of units. The accelerogram is scaled 
uniformly to achieve the specified peak values of GMAXH and GMAXV. Since data is 
read in free format, as many lines as necessary to read the entire wave must be input. 


GO TO DATA SET U 


e USER_TEXT Reference information: upto 80 characters of text 
e ICNTRL Cyclic Analysis option 
= 0, Force controlled input 


= 1, displacement controlled input 


e NLDED Number of story levels at which the force /displacement is applied 
e NSTLD(),I=1,NLDED List of story levels at which the force or displacement is applied 

e NPTS Number of points to be read in force or displacement history 

e F(I,1),I=1,NPTS first data set (NPTS) at story level NSTLD(1) 
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e F(I,2),I=1,NPTS next data set (NPTS) at story level NSTLD(2) 
e (repeat for each of NLDED levels) 
e ITCAL No. of points to interpolate between prescribed load steps 


The analysis is performed at ITCAL interpolated points for each step 


e USER_TEXT Reference information: upto 80 characters of text 
e NSOUT,DTOUT,ISO(),I=1,NSOUT 
e FNAME (1) 


e (continue with filenames for each of NSOUT output sets) 


NSOUT = No of output histories 

DTOUT = Output time interval 

ISO(I) = Output story numbers 

FNAME(I) = Filename to store time history output for story number ISO(1) 


NOTES: For the quasi-static cyclic analysis option, DTOUT refers to the number of steps 
between output printing; for example, DTOUT=2 will print results every 2 steps. 


e USER_TEXT Reference information: upto 80 characters 
e KCOUT, KBOUT, KWOUT, KSOUT, KPOUT 


KCOUT = Number of columns for which hysteresis output is required 
KBOUT = Number of beams for which hysteresis output is required 
KWOUT = Number of walls for which hysteresis output is required 
KSOUT = Number of springs for which hysteresis output is required 
KPOUT = Number of panels for which hysteresis output is required 


COLUMN OUTPUT SPECIFICATION 

SKIP THIS INPUT IF KCOUT = 0 

e USER_TEXT Reference information: upto 80 characters 

e ICLIST(), I=1,KCOUT List of column numbers for which 
moment-curvature hysteresis is required 


BEAM OUTPUT SPECIFICATION 

SKIP THIS INPUT IF KBOUT = 0 

e USER_TEXT Reference information: upto 80 characters 

e IBLIST(D, I=1,KBOUT List of beam numbers for which 
moment-curvature hysteresis is required 
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SHEAR WALL OUTPUT SPECIFICATION 
SKIP THIS INPUT IF KWOUT = 0 


e USER_TEXT Reference information: upto 80 characters 

e IWLIST(D, I=1,KWOUT List of shear wall numbers for which 
moment-curvature and shear-strain hysteresis 
is required 


DISCRETE SPRING OUTPUT SPECIFICATION 

SKIP THIS INPUT IF KSOUT = 0 

e USER_TEXT Reference information: upto 80 characters 

e ISLIST(), I=1,KSOUT List of spring numbers for which 
moment-rotation hysteresis is required 


JOINT PANEL OUTPUT SPECIFICATION 
SKIP THIS INPUT IF KPOUT = 0 
e USER_TEXT Reference information: upto 80 characters 
e ISLIST(), I=1,KPOUT List of joint panel numbers for which 
shear vs. panel deformation hysteresis is required 


NOTES: All the output generated in this section refers to moment-curvature hysteresis for 
beams, columns and shear-walls; in addition shear vs. shear strain history is generated 
for walls; whereas moment-rotation hysteresis is produced for the discrete spring 
elements. Output filenames are generated as follows: 

IF KCOUT = 2, AND ICLIST(1) = 3 AND ICLIST(2) = 12, THEN THE FOLLOWING 
FILES WILL BE CREATED: 

COL_003.PRN and COL_012.PRN 

(where 3 and 12 refer to the element numbers for which output is requested) 
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APPENDIX B 


SAMPLE DATA SETS 
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DATA SET 


SAMPLE PROBLEM TO VERIFY MEMBER MODEL : NONLINEAR STATIC ANALYSIS 
CONTROL DATA 
hy Oe Ep Es a 
ELEMENT TYPES 
2100000 
ELEMENT DATA 
21000000 
UNIT SYSTEM 
a | 
FLOOR ELEVATIONS 
120.0 
DESCRIPTION OF IDENTICAL FRAMES 
1 
PLAN CONFIGURATION 
2 
NODAL WEIGHTS 
11 100.0 100.0 
ENVELOPE GENERATION 
1 
HYSTERESIS MODELLING 
1 
1 2 0.05 0.0 0.0 0.0 0.0 0.0 
COLUMN PROPERTIES 
4 
MOMENT CURVATURE ENVELOPE FOR THE STEEL COLUMN 
1 120.0 0.0 0.0 
1 7.5E6 5000.0 0.0 1800.0 1800.0 
1 7.5E6 5000.0 0.0 1800.0 1800.0 
2 120.0 0.0 0.0 
1 7.5E6 5000.0 0.0 1800.0 1800.0 
1 7.5E6 5000.0 0.0 1800.0 1800.0 
BEAM PROPERTIES 
4 
MOMENT CURVATURE ENVELOPE FOR BEAM 
1 180.0 0.0 0.0 
1 5E6 0.0 650.0 650.0 
1 5E6 0.0 650.0 650.0 
COLUMN CONNECTIONS 
2182. 1,031 
feanlae2e0 1 
BEAM CONNECTIONS 
22191 2132 
ANALYSIS TYPE INELASTIC INCREMENTAL ANALYSIS WITH STATIC LOADS 
1 
STATIC ANALYSIS - LATERAL LOAD AT FLOOR 
0100 
20 1 
LATERALLY LOADED JOINTS (Change magnitude for other cases) 
Lelel 4530 


SAMPLE PROBLEM TO VERIFY MEMBER MODEL : NONLINEAR DYNAMIC ANALYSIS 


CONTROL DATA 

Wega ho he a 

ELEMENT TYPES 

anie0. (0, 07020 

ELEMENT DATA 

2517070, 0 070.0 

UNIT SYSTEM 

1 

FLOOR ELEVATIONS 

120.0 

DESCRIPTION OF IDENTICAL FRAMES 

1 

PLAN CONFIGURATION 

2 

NODAL WEIGHTS 

11 100.0 100.0 

ENVELOPE GENERATION 

1 

HYSTERESIS MODELLING 

ak 

12 0.05 0.0 0.0 0.0 0.0 0.0 

COLUMN PROPERTIES 

4 

MOMENT CURVATURE ENVELOPE FOR THE STEEL COLUMN 

1 120.0 0.0 0.0 
1 7.5E6 5000.0 0.0 1800.0 1800.0 
1 7.5E6 5000.0 0.0 1800.0 1800.0 

2 120.0 0.0 0.0 
1 7.5E6 5000.0 0.0 1800.0 1800.0 
1 7.5E6 5000.0 0.0 1800.0 1800.0 

BEAM PROPERTIES 

4 

MOMENT CURVATURE ENVELOPE FOR BEAM 

1 180.0 0.0 0.0 

1 5E6 0.0 650.0 650.0 

1 5E6 0.0 650.0 650.0 

COLUMN CONNECTIONS 

pbk aR ah. to) at 

Zea. 42021 

BEAM CONNECTIONS 

Lie ed ek ee 

ANALYSIS TYPE INELASTIC INCREMENTAL ANALYSIS WITH STATIC LOADS 

3 

Long term loads 

0000 

Dynamic analysis 

0.12 0.0 0.02 10.0 0.0 

Wave Data 

0, 1001, 0.02 

ElCentro 

ELC .DAT 

OUTPUT CONTROL 

0 204m) 

ID4G. PRN 

ELEMENT HYSTERESIS OUTPUT 

2. 1-030 

COLUMN OUTPUT 

Lae 

BEAM OUTPUT 

oo 
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Lehigh Test - Panel Zone Deformations 
Control data 
Behe, Oe ag 0 


Elem types 

eee 0 pe Ot Os a 
Elem data 

ay ease 0s 0, OF LG 0 
Units 

1 

Floor elev 

54.0 108.0 
Duplicate frames 

ui 

Plan config 

3 


Nodal weights 

1,, 4, 10,0 10.0°10.0 

Qf the? -O.Onr10 70, 0<0 

Env generation 

0 

Steel prop 

1, 40.0, 58.8, 30000.0, 300.0, 3.0 
Hys model 

2 

1, 2, 0.01 00000 
Column prop 

3 

Steel section data 
Veetk tiyaeoe. 0. 0.0712: 
2, 1, 1, S@30 1220. C. 
Beam prop 

3 

Steel section dat 
4 ods 08.0; 0% 
2, 1, 1, 68.@0, <7? 
Joint Panels 

11 24.0 14.0 0.44 11000.0 39.0 39.0 
Column conn 

LiL dy a2 07 1 

af ares as 17 2 

Beam conn 

apes hs by ky 

hy a? ee i ae 

Panel Location 


-16 999.0 357 360 
-16 999.0 357 360 


om ©. 
uu 
nan 


0 18.2 10.21 1550.0 131.0 153.0 


a 
07.0000 
00.0 000 0 18.2 10.21 1550.0 131.0 153.0 


a Roe 1a ee amb 

Analysis 

4 

Long Term Loads - none 

O,. 0, 0,.0 

Quasistatic Loading : Disp Control 

1 

2 

27 #! Refers to DOFs 

32 

0 -2 -4 -6 -6.65 -6 -4 -2 0 2 4 6 6.70.6 C45 2.20 
-2 -4 -6 -7.75 -6 -4 -2 0 2 4 6 10.2 6 4 2 

0 2 4 6 6.65 6 4 2 0 -2 -4 -6 -6.7 -6 -4 -2 0O 
FE WS Bo f=) 6 4 2 O -2 -4 -6 -10.2 -6 -4 -2 

40 

Output control 

ively a 


LEHI. PRN Next 2 lines: Misc Output - None requested ; 0000 0 


B-3 


Verification of Hysteresis Model 
CONTROL DATA 


UCB TEST : 


2, 1, 90, 2, 
ELEM TYPES 
2, 1, 0, O, 
ELEM DATA 
3, 2, 0, 
UNITS 

ak 

FLOOR ELEV 

17.0, 158.0 
DUPLICATE FRAMES 

a 

PLAN CONFIG 

3 

NODAL WEIGHTS 

1, 2, 210.0°10.0 10.0 
2, 1, 0.0 10.0 0.0 
ENV GENERATION 
0 

STEEL PROP 

1 58.5,..58-8, 
2753.57.02 .8, 
HYS MODEL 

1 

1, 3, 0.2 0.4 0.24 0.4 0.9 1.2 
COLUMN PROP 

3 

STEEL SECTION DATA 

i, 1, 1, 17.090 .0°0. 
2, 1, 1, 141.0 0.0 O. 
BEAM PROP 

3 

STEEL SECTION DATA 

1, 1, 1, 68.0, 
COLUMN CONN 
L, @ige Le; 
a5 wee 17s, 
CF a re ay 
BEAM CONN 
ah a Bees 
ey Ls 
ANALYSIS 
4 

STATIC 
O75.0, 70, 0 
QUASISTATIC 


0 
0, 0, O 


0, O, 


30000.0, 
30000.0, 


5.0, 


-0 
5.0, 2.0 


0 0.0 0000 0 
0 0.0 0000 


1, 
1, 


wre NEB 


FOOOCOUrRNE 


OUTPUT CONTROL 
Licgol ,ég@ 

BEAM_TIP.PRN 
MISC OUTPUT; 


Next Line: 0 0000 


B-4 


0 


30.6 100.0 2000 3000 
44.2 22.4 9040.0 393.0 581.0 


0.0 0.0 0 00 0 75.6 19.25 3400.0 415.0 487.0 


“1.5 ~-1.75 
-1.0 0.0 
1.85. 125.1750 


ii, | 


Bru 
mM 


a 
as 2 
ae - 


oy. 


, 


